3.3.17 \(\int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx\) [217]

3.3.17.1 Optimal result
3.3.17.2 Mathematica [A] (verified)
3.3.17.3 Rubi [A] (verified)
3.3.17.4 Maple [B] (verified)
3.3.17.5 Fricas [F]
3.3.17.6 Sympy [F]
3.3.17.7 Maxima [F]
3.3.17.8 Giac [F(-2)]
3.3.17.9 Mupad [F(-1)]

3.3.17.1 Optimal result

Integrand size = 29, antiderivative size = 314 \[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=-\frac {b^2 c^2 \sqrt {d-c^2 d x^2}}{3 x}-\frac {b^2 c^3 \sqrt {d-c^2 d x^2} \arcsin (c x)}{3 \sqrt {1-c^2 x^2}}-\frac {b c \sqrt {1-c^2 x^2} \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))}{3 x^2}+\frac {i c^3 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}-\frac {2 b c^3 \sqrt {d-c^2 d x^2} (a+b \arcsin (c x)) \log \left (1-e^{2 i \arcsin (c x)}\right )}{3 \sqrt {1-c^2 x^2}}+\frac {i b^2 c^3 \sqrt {d-c^2 d x^2} \operatorname {PolyLog}\left (2,e^{2 i \arcsin (c x)}\right )}{3 \sqrt {1-c^2 x^2}} \]

output
-1/3*(-c^2*d*x^2+d)^(3/2)*(a+b*arcsin(c*x))^2/d/x^3-1/3*b^2*c^2*(-c^2*d*x^ 
2+d)^(1/2)/x-1/3*b^2*c^3*arcsin(c*x)*(-c^2*d*x^2+d)^(1/2)/(-c^2*x^2+1)^(1/ 
2)+1/3*I*c^3*(a+b*arcsin(c*x))^2*(-c^2*d*x^2+d)^(1/2)/(-c^2*x^2+1)^(1/2)-2 
/3*b*c^3*(a+b*arcsin(c*x))*ln(1-(I*c*x+(-c^2*x^2+1)^(1/2))^2)*(-c^2*d*x^2+ 
d)^(1/2)/(-c^2*x^2+1)^(1/2)+1/3*I*b^2*c^3*polylog(2,(I*c*x+(-c^2*x^2+1)^(1 
/2))^2)*(-c^2*d*x^2+d)^(1/2)/(-c^2*x^2+1)^(1/2)-1/3*b*c*(a+b*arcsin(c*x))* 
(-c^2*x^2+1)^(1/2)*(-c^2*d*x^2+d)^(1/2)/x^2
 
3.3.17.2 Mathematica [A] (verified)

Time = 1.18 (sec) , antiderivative size = 248, normalized size of antiderivative = 0.79 \[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=\frac {\sqrt {d-c^2 d x^2} \left (2 b^2 \left (i c^3 x^3-\sqrt {1-c^2 x^2}+c^2 x^2 \sqrt {1-c^2 x^2}\right ) \arcsin (c x)^2-b \arcsin (c x) \left (2 b c x+3 a \sqrt {1-c^2 x^2}+a \cos (3 \arcsin (c x))+4 b c^3 x^3 \log \left (1-e^{2 i \arcsin (c x)}\right )\right )-2 \left (a b c x+b^2 c^2 x^2 \sqrt {1-c^2 x^2}+a^2 \left (1-c^2 x^2\right )^{3/2}+2 a b c^3 x^3 \log (c x)\right )+2 i b^2 c^3 x^3 \operatorname {PolyLog}\left (2,e^{2 i \arcsin (c x)}\right )\right )}{6 x^3 \sqrt {1-c^2 x^2}} \]

input
Integrate[(Sqrt[d - c^2*d*x^2]*(a + b*ArcSin[c*x])^2)/x^4,x]
 
output
(Sqrt[d - c^2*d*x^2]*(2*b^2*(I*c^3*x^3 - Sqrt[1 - c^2*x^2] + c^2*x^2*Sqrt[ 
1 - c^2*x^2])*ArcSin[c*x]^2 - b*ArcSin[c*x]*(2*b*c*x + 3*a*Sqrt[1 - c^2*x^ 
2] + a*Cos[3*ArcSin[c*x]] + 4*b*c^3*x^3*Log[1 - E^((2*I)*ArcSin[c*x])]) - 
2*(a*b*c*x + b^2*c^2*x^2*Sqrt[1 - c^2*x^2] + a^2*(1 - c^2*x^2)^(3/2) + 2*a 
*b*c^3*x^3*Log[c*x]) + (2*I)*b^2*c^3*x^3*PolyLog[2, E^((2*I)*ArcSin[c*x])] 
))/(6*x^3*Sqrt[1 - c^2*x^2])
 
3.3.17.3 Rubi [A] (verified)

Time = 0.82 (sec) , antiderivative size = 206, normalized size of antiderivative = 0.66, number of steps used = 13, number of rules used = 12, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.414, Rules used = {5186, 5190, 247, 223, 5136, 3042, 25, 4200, 25, 2620, 2715, 2838}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx\)

\(\Big \downarrow \) 5186

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \int \frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{x^3}dx}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 5190

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \left (c^2 \left (-\int \frac {a+b \arcsin (c x)}{x}dx\right )+\frac {1}{2} b c \int \frac {\sqrt {1-c^2 x^2}}{x^2}dx-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}\right )}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 247

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \left (c^2 \left (-\int \frac {a+b \arcsin (c x)}{x}dx\right )+\frac {1}{2} b c \left (c^2 \left (-\int \frac {1}{\sqrt {1-c^2 x^2}}dx\right )-\frac {\sqrt {1-c^2 x^2}}{x}\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}\right )}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 223

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \left (c^2 \left (-\int \frac {a+b \arcsin (c x)}{x}dx\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 5136

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \left (c^2 \left (-\int \frac {\sqrt {1-c^2 x^2} (a+b \arcsin (c x))}{c x}d\arcsin (c x)\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \left (c^2 \left (-\int -\left ((a+b \arcsin (c x)) \tan \left (\arcsin (c x)+\frac {\pi }{2}\right )\right )d\arcsin (c x)\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 25

\(\displaystyle \frac {2 b c \sqrt {d-c^2 d x^2} \left (c^2 \int (a+b \arcsin (c x)) \tan \left (\arcsin (c x)+\frac {\pi }{2}\right )d\arcsin (c x)-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}-\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}\)

\(\Big \downarrow \) 4200

\(\displaystyle -\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}+\frac {2 b c \sqrt {d-c^2 d x^2} \left (-\left (c^2 \left (2 i \int -\frac {e^{2 i \arcsin (c x)} (a+b \arcsin (c x))}{1-e^{2 i \arcsin (c x)}}d\arcsin (c x)-\frac {i (a+b \arcsin (c x))^2}{2 b}\right )\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}\)

\(\Big \downarrow \) 25

\(\displaystyle -\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}+\frac {2 b c \sqrt {d-c^2 d x^2} \left (-\left (c^2 \left (-2 i \int \frac {e^{2 i \arcsin (c x)} (a+b \arcsin (c x))}{1-e^{2 i \arcsin (c x)}}d\arcsin (c x)-\frac {i (a+b \arcsin (c x))^2}{2 b}\right )\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}\)

\(\Big \downarrow \) 2620

\(\displaystyle -\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}+\frac {2 b c \sqrt {d-c^2 d x^2} \left (-\left (c^2 \left (-2 i \left (\frac {1}{2} i \log \left (1-e^{2 i \arcsin (c x)}\right ) (a+b \arcsin (c x))-\frac {1}{2} i b \int \log \left (1-e^{2 i \arcsin (c x)}\right )d\arcsin (c x)\right )-\frac {i (a+b \arcsin (c x))^2}{2 b}\right )\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}\)

\(\Big \downarrow \) 2715

\(\displaystyle -\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}+\frac {2 b c \sqrt {d-c^2 d x^2} \left (-\left (c^2 \left (-2 i \left (\frac {1}{2} i \log \left (1-e^{2 i \arcsin (c x)}\right ) (a+b \arcsin (c x))-\frac {1}{4} b \int e^{-2 i \arcsin (c x)} \log \left (1-e^{2 i \arcsin (c x)}\right )de^{2 i \arcsin (c x)}\right )-\frac {i (a+b \arcsin (c x))^2}{2 b}\right )\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}\)

\(\Big \downarrow \) 2838

\(\displaystyle -\frac {\left (d-c^2 d x^2\right )^{3/2} (a+b \arcsin (c x))^2}{3 d x^3}+\frac {2 b c \sqrt {d-c^2 d x^2} \left (-\left (c^2 \left (-2 i \left (\frac {1}{2} i \log \left (1-e^{2 i \arcsin (c x)}\right ) (a+b \arcsin (c x))+\frac {1}{4} b \operatorname {PolyLog}\left (2,e^{2 i \arcsin (c x)}\right )\right )-\frac {i (a+b \arcsin (c x))^2}{2 b}\right )\right )-\frac {\left (1-c^2 x^2\right ) (a+b \arcsin (c x))}{2 x^2}+\frac {1}{2} b c \left (-c \arcsin (c x)-\frac {\sqrt {1-c^2 x^2}}{x}\right )\right )}{3 \sqrt {1-c^2 x^2}}\)

input
Int[(Sqrt[d - c^2*d*x^2]*(a + b*ArcSin[c*x])^2)/x^4,x]
 
output
-1/3*((d - c^2*d*x^2)^(3/2)*(a + b*ArcSin[c*x])^2)/(d*x^3) + (2*b*c*Sqrt[d 
 - c^2*d*x^2]*(-1/2*((1 - c^2*x^2)*(a + b*ArcSin[c*x]))/x^2 + (b*c*(-(Sqrt 
[1 - c^2*x^2]/x) - c*ArcSin[c*x]))/2 - c^2*(((-1/2*I)*(a + b*ArcSin[c*x])^ 
2)/b - (2*I)*((I/2)*(a + b*ArcSin[c*x])*Log[1 - E^((2*I)*ArcSin[c*x])] + ( 
b*PolyLog[2, E^((2*I)*ArcSin[c*x])])/4))))/(3*Sqrt[1 - c^2*x^2])
 

3.3.17.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 223
Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Simp[ArcSin[Rt[-b, 2]*(x/Sqrt 
[a])]/Rt[-b, 2], x] /; FreeQ[{a, b}, x] && GtQ[a, 0] && NegQ[b]
 

rule 247
Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(c*x)^ 
(m + 1)*((a + b*x^2)^p/(c*(m + 1))), x] - Simp[2*b*(p/(c^2*(m + 1)))   Int[ 
(c*x)^(m + 2)*(a + b*x^2)^(p - 1), x], x] /; FreeQ[{a, b, c}, x] && GtQ[p, 
0] && LtQ[m, -1] &&  !ILtQ[(m + 2*p + 3)/2, 0] && IntBinomialQ[a, b, c, 2, 
m, p, x]
 

rule 2620
Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/ 
((a_) + (b_.)*((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)), x_Symbol] :> Simp 
[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x] - Si 
mp[d*(m/(b*f*g*n*Log[F]))   Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x 
)))^n/a)], x], x] /; FreeQ[{F, a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]
 

rule 2715
Int[Log[(a_) + (b_.)*((F_)^((e_.)*((c_.) + (d_.)*(x_))))^(n_.)], x_Symbol] 
:> Simp[1/(d*e*n*Log[F])   Subst[Int[Log[a + b*x]/x, x], x, (F^(e*(c + d*x) 
))^n], x] /; FreeQ[{F, a, b, c, d, e, n}, x] && GtQ[a, 0]
 

rule 2838
Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2 
, (-c)*e*x^n]/n, x] /; FreeQ[{c, d, e, n}, x] && EqQ[c*d, 1]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4200
Int[((c_.) + (d_.)*(x_))^(m_.)*tan[(e_.) + Pi*(k_.) + (f_.)*(x_)], x_Symbol 
] :> Simp[I*((c + d*x)^(m + 1)/(d*(m + 1))), x] - Simp[2*I   Int[(c + d*x)^ 
m*E^(2*I*k*Pi)*(E^(2*I*(e + f*x))/(1 + E^(2*I*k*Pi)*E^(2*I*(e + f*x)))), x] 
, x] /; FreeQ[{c, d, e, f}, x] && IntegerQ[4*k] && IGtQ[m, 0]
 

rule 5136
Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_.)/(x_), x_Symbol] :> Subst[Int[( 
a + b*x)^n*Cot[x], x], x, ArcSin[c*x]] /; FreeQ[{a, b, c}, x] && IGtQ[n, 0]
 

rule 5186
Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))^(n_.)*((f_.)*(x_))^(m_)*((d_) + (e_. 
)*(x_)^2)^(p_), x_Symbol] :> Simp[(f*x)^(m + 1)*(d + e*x^2)^(p + 1)*((a + b 
*ArcSin[c*x])^n/(d*f*(m + 1))), x] - Simp[b*c*(n/(f*(m + 1)))*Simp[(d + e*x 
^2)^p/(1 - c^2*x^2)^p]   Int[(f*x)^(m + 1)*(1 - c^2*x^2)^(p + 1/2)*(a + b*A 
rcSin[c*x])^(n - 1), x], x] /; FreeQ[{a, b, c, d, e, f, m, p}, x] && EqQ[c^ 
2*d + e, 0] && GtQ[n, 0] && EqQ[m + 2*p + 3, 0] && NeQ[m, -1]
 

rule 5190
Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))*((f_.)*(x_))^(m_)*((d_) + (e_.)*(x_) 
^2)^(p_.), x_Symbol] :> Simp[(f*x)^(m + 1)*(d + e*x^2)^p*((a + b*ArcSin[c*x 
])/(f*(m + 1))), x] + (-Simp[b*c*(d^p/(f*(m + 1)))   Int[(f*x)^(m + 1)*(1 - 
 c^2*x^2)^(p - 1/2), x], x] - Simp[2*e*(p/(f^2*(m + 1)))   Int[(f*x)^(m + 2 
)*(d + e*x^2)^(p - 1)*(a + b*ArcSin[c*x]), x], x]) /; FreeQ[{a, b, c, d, e, 
 f}, x] && EqQ[c^2*d + e, 0] && IGtQ[p, 0] && ILtQ[(m + 1)/2, 0]
 
3.3.17.4 Maple [B] (verified)

Both result and optimal contain complex but leaf count of result is larger than twice the leaf count of optimal. 2040 vs. \(2 (294 ) = 588\).

Time = 0.29 (sec) , antiderivative size = 2041, normalized size of antiderivative = 6.50

method result size
default \(\text {Expression too large to display}\) \(2041\)
parts \(\text {Expression too large to display}\) \(2041\)

input
int((-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2/x^4,x,method=_RETURNVERBOSE)
 
output
-1/3*a^2/d/x^3*(-c^2*d*x^2+d)^(3/2)-1/3*I*b^2*(-d*(c^2*x^2-1))^(1/2)/(3*c^ 
4*x^4-3*c^2*x^2+1)/(c^2*x^2-1)*(-c^2*x^2+1)^(1/2)*c^3-2*I*b^2*(-d*(c^2*x^2 
-1))^(1/2)*(-c^2*x^2+1)^(1/2)*c^3/(3*c^2*x^2-3)*arcsin(c*x)^2-2*I*b^2*(-d* 
(c^2*x^2-1))^(1/2)*(-c^2*x^2+1)^(1/2)*c^3/(3*c^2*x^2-3)*polylog(2,-I*c*x-( 
-c^2*x^2+1)^(1/2))-2*I*b^2*(-d*(c^2*x^2-1))^(1/2)*(-c^2*x^2+1)^(1/2)*c^3/( 
3*c^2*x^2-3)*polylog(2,I*c*x+(-c^2*x^2+1)^(1/2))+2*b^2*(-d*(c^2*x^2-1))^(1 
/2)*(-c^2*x^2+1)^(1/2)*c^3/(3*c^2*x^2-3)*arcsin(c*x)*ln(1+I*c*x+(-c^2*x^2+ 
1)^(1/2))+2*b^2*(-d*(c^2*x^2-1))^(1/2)*(-c^2*x^2+1)^(1/2)*c^3/(3*c^2*x^2-3 
)*arcsin(c*x)*ln(1-I*c*x-(-c^2*x^2+1)^(1/2))+1/3*b^2*(-d*(c^2*x^2-1))^(1/2 
)/(3*c^4*x^4-3*c^2*x^2+1)*x^3/(c^2*x^2-1)*(-c^2*x^2+1)*c^6-b^2*(-d*(c^2*x^ 
2-1))^(1/2)/(3*c^4*x^4-3*c^2*x^2+1)/(c^2*x^2-1)*arcsin(c*x)*(-c^2*x^2+1)^( 
1/2)*c^3+b^2*(-d*(c^2*x^2-1))^(1/2)/(3*c^4*x^4-3*c^2*x^2+1)*x^5/(c^2*x^2-1 
)*arcsin(c*x)^2*c^8-3*b^2*(-d*(c^2*x^2-1))^(1/2)/(3*c^4*x^4-3*c^2*x^2+1)*x 
^3/(c^2*x^2-1)*arcsin(c*x)^2*c^6+10/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(3*c^4*x^ 
4-3*c^2*x^2+1)*x/(c^2*x^2-1)*arcsin(c*x)^2*c^4-5/3*b^2*(-d*(c^2*x^2-1))^(1 
/2)/(3*c^4*x^4-3*c^2*x^2+1)/x/(c^2*x^2-1)*arcsin(c*x)^2*c^2+I*b^2*(-d*(c^2 
*x^2-1))^(1/2)/(3*c^4*x^4-3*c^2*x^2+1)*x^4/(c^2*x^2-1)*(-c^2*x^2+1)^(1/2)* 
arcsin(c*x)^2*c^7-1/3*I*b^2*(-d*(c^2*x^2-1))^(1/2)/(3*c^4*x^4-3*c^2*x^2+1) 
*x^3/(c^2*x^2-1)*(-c^2*x^2+1)*arcsin(c*x)*c^6-I*b^2*(-d*(c^2*x^2-1))^(1/2) 
/(3*c^4*x^4-3*c^2*x^2+1)*x^2/(c^2*x^2-1)*(-c^2*x^2+1)^(1/2)*arcsin(c*x)...
 
3.3.17.5 Fricas [F]

\[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=\int { \frac {\sqrt {-c^{2} d x^{2} + d} {\left (b \arcsin \left (c x\right ) + a\right )}^{2}}{x^{4}} \,d x } \]

input
integrate((-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2/x^4,x, algorithm="frica 
s")
 
output
integral(sqrt(-c^2*d*x^2 + d)*(b^2*arcsin(c*x)^2 + 2*a*b*arcsin(c*x) + a^2 
)/x^4, x)
 
3.3.17.6 Sympy [F]

\[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=\int \frac {\sqrt {- d \left (c x - 1\right ) \left (c x + 1\right )} \left (a + b \operatorname {asin}{\left (c x \right )}\right )^{2}}{x^{4}}\, dx \]

input
integrate((-c**2*d*x**2+d)**(1/2)*(a+b*asin(c*x))**2/x**4,x)
 
output
Integral(sqrt(-d*(c*x - 1)*(c*x + 1))*(a + b*asin(c*x))**2/x**4, x)
 
3.3.17.7 Maxima [F]

\[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=\int { \frac {\sqrt {-c^{2} d x^{2} + d} {\left (b \arcsin \left (c x\right ) + a\right )}^{2}}{x^{4}} \,d x } \]

input
integrate((-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2/x^4,x, algorithm="maxim 
a")
 
output
1/3*((-1)^(-2*c^2*d*x^2 + 2*d)*c^2*d^(3/2)*log(-2*c^2*d + 2*d/x^2) + c^2*d 
^(3/2)*log(x^2 - 1/c^2) - sqrt(c^4*d*x^4 - 2*c^2*d*x^2 + d)*d/x^2)*a*b*c/d 
 - 2/3*(-c^2*d*x^2 + d)^(3/2)*a*b*arcsin(c*x)/(d*x^3) + 1/3*((c^2*x^2 - 1) 
*sqrt(c*x + 1)*sqrt(-c*x + 1)*sqrt(d)*arctan2(c*x, sqrt(c*x + 1)*sqrt(-c*x 
 + 1))^2 - sqrt(d)*x^3*integrate(2*(c^3*x^2 - c)*arctan2(c*x, sqrt(c*x + 1 
)*sqrt(-c*x + 1))/x^3, x))*b^2/x^3 - 1/3*(-c^2*d*x^2 + d)^(3/2)*a^2/(d*x^3 
)
 
3.3.17.8 Giac [F(-2)]

Exception generated. \[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=\text {Exception raised: TypeError} \]

input
integrate((-c^2*d*x^2+d)^(1/2)*(a+b*arcsin(c*x))^2/x^4,x, algorithm="giac" 
)
 
output
Exception raised: TypeError >> an error occurred running a Giac command:IN 
PUT:sage2:=int(sage0,sageVARx):;OUTPUT:sym2poly/r2sym(const gen & e,const 
index_m & i,const vecteur & l) Error: Bad Argument Value
 
3.3.17.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {d-c^2 d x^2} (a+b \arcsin (c x))^2}{x^4} \, dx=\int \frac {{\left (a+b\,\mathrm {asin}\left (c\,x\right )\right )}^2\,\sqrt {d-c^2\,d\,x^2}}{x^4} \,d x \]

input
int(((a + b*asin(c*x))^2*(d - c^2*d*x^2)^(1/2))/x^4,x)
 
output
int(((a + b*asin(c*x))^2*(d - c^2*d*x^2)^(1/2))/x^4, x)